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We have used single-dish centimeter-band, multi-frequency linear polarization and total flux density monitoring 
observations obtained by UMRAO at 14.5, 8, and 4.8 GHz during 7-ray flaring to constrain propagating shock- 
in-jet models; this procedure allows us to probe jet conditions at or near the presumed parsec-scale 7-ray emission 
site within the relativistic jet. Results are presented for the blazar 0420-014 during the 7-ray flare detected by 
the LAT which peaked in late January 2010. Our work identifies the passage of three forward-moving transverse 
shocks during the radio-band flaring, a shock Lorentz factor of 12, and an observer's viewing angle of 2° and 
sets limits on the energy distribution of the radiating particles. 



1 . Overview 

To investigate whether shocks are present during 
7-ray flares, and to derive source parameters during 
7-ray flaring we have monitored ~30 blazars for to- 
tal flux density and linear polarization (hereafter LP) 
with the University of Michigan 26-m telescope (UM- 
RAO) at 14.5, 8.0, and 4.8 GHz. Specific goals of 
our work are to determine whether evidence for the 
occurrence of internal shocks can be found in the LP 
data; if present, we use comparisons of the UMRAO 
data with simulations from radiative transfer model- 
ing incorporating propagating shocks to constrain jet 
flow conditions during the 7-ray flaring. We present 
here results for the HPQ 0420-014, analyzed as part 
of a pilot study. This source was selected on the basis 
of showing both strong and resolved radio band flares 
since 2008.5, and 7-ray flaring at GeV energies in the 
fermi-LAT monitoring data. 



2. Historical Character of the Radio-band 
Variability and Relevance to 7-ray Flaring 

Centimeter-band total flux density (S) and lin- 
ear polarization (LP) monitoring observations of 
the 7-ray-bright blazar 0420-014 (PKS 0420-01; 
2FGLJ0423.2-0120) have been obtained with the Uni- 
versity of Michigan (UMRAO) 26-m radio telescope 
since the late 1970s. These data, as shown in Fig- 
ure [TJ document that the centimeter-band variabil- 
ity has been intense and continuous during the past 
3.5 decades. A peak amplitude of Swl4 Jy at 14.5 
GHz was reached in late 2003, making this source 
one of the brightest blazars historically in the cen- 
timeter band. Prior to this large outburst, the 
spectrum was very flat, but the spectrum steep- 
ened during the subsequent large flare. Temporally- 
resolved flares in fractional LP have also occurred 



throughout the time period monitored; these ob- 
servations identify a peak amplitude of 5%. The 
electric vector position angle (EVPA: shown in the 
top panel) has exhibited a series of sharp, system- 
atic changes lasting from several months to years. 
Complementary, long-term spatially-resolved VLBI 
imaging data from the MO J AVE program (at 15.4 
GHz) have been obtained since the mid 1990s with 
a typical cadence of a few observations per year; 
these imaging data track the motion of 3 compo- 
nents from 1995 through 2008 and identify a maxi- 
mum apparent component speed of /3app=7.36 c (see 
http:www.physics.purdue.edu/astro/MOJAVE). In a 
map obtained in February 2010, a new, very bright 
component is apparent; this feature is moving at a 
different position angle from its predecessors and with 
/3app=3.99±0.80 c based on 6 epochs of observation [l[. 

The source was detected by EGRET in the 1990s, 
with a highly significant detection (TS=46) in the 
viewing period centered on 1992.16 during a 
plateau in a centimeter-band fiare. Shortly after 
the launch of Fermi the source flared in the radio- 
band, and contemporaneous 7-ray-radio-band flaring 
was detected by the LAT. 

An explanation for the cause of the radio-band flar- 
ing commonly accepted since the mid-1980s is internal 
shocks in the parsec-scale region of the jet where the 
synchrotron emission arises Q, [4|. These structures 
have been shown to develop naturally in hydrody- 
namic simulations of AGN jets Q, validating this sce- 
nario. Several radio band studies based on the devel- 
opment of correlated cross-band activity support the 
idea that the 7-ray emission is also produced within 
the parsec-scale region of the jet, at or near the 43 
GHz 'core', e.g. Q. Thus, the radio band data probe 
the region in which the 7-ray flares are plausibly pro- 
duced. Evidence for the presence of shocks is of inter- 
est since shocks are a potential mechanism for acceler- 
ating the emitting particles producing the synchrotron 
emission to the high energies required for the produc- 
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Figure 1: Two- week averages of the long-term UMRAO 
data showing from bottom to top: S, fractional LP (P%), 
and EVPA. The observations at the three frequencies are 
symbol and color coded as indicated in the upper left. 
The source-integrated rotation measure of -13 rad/m^ 
determined from 18-20.5 cm VLA data 8] produces a 
negligible shift in the UMRAO EVPAs. " 



tion of 7-ray flares. An earlier study of our target by 
Stevens et al. [7]) compared the spectral evolution 
of total flux density flares during the period 1989.0- 
1994, which includes the strong EGRET detection, 
with the predictions of the Marscher & Gear Q shock 
model using high frequency data (375 - 22 GHz); the 
spectral behavior was found to be consistent with the 
shock model predictions. However, discrepancies be- 
tween the observations and the expected spectral evo- 
lution were also identified. These were qualitatively 
attributed to bending in an initially-straight jet ori- 
ented close to the line of sight. 



3. Description of The Model 

In our model, we envisage a propagating shock pass- 
ing through a medium containing a turbulent mag- 
netic field. As a consequence of the compression as- 
sociated with the passage of the shock, there is an 
increase in the fractional linear polarization, and an 
increase in the particle density (and, hence, the emis- 
sivity). The outburst is associated with the prop- 
agation of the region bounded by the limits of the 
shocked flow. In an earlier phase of our program and 
following the method of ^ we developed propagating 
shock models allowing for the shock to be oriented 
at an arbitrary angle to the jet flow direction. This 
work is described in [l^]. Each shock is specified by 
its sense (forward or reverse), fractional length rela- 
tive to the quiescent flow, compression factor, and the 
shock's obliquity (the orientation relative to the flow 
direction) . The choice of value for an additional free 
parameter specifying the observer's orientation rela- 
tive to the shock plane was shown to have little af- 



+ ' t +' ' 

t • J + 
. 1 "t 

\ «t r 



' ii + ■ 

1 , , t 


(ill 1 ' 


0420-01 4 
^ ' ' J i ' 




s « ^ ^ i 




1 1 1 1 1 1 





5000 5200 5400 5600 

Julian Date - 2450000 



Figure 2: The 7-ray and centimeter band light curves 
during the radio-band event modeled. The time window 
shown is March 9, 2009 - February 7, 2011. From 
bottom to top: (panel 1) the weekly-binned 7-ray light 
curve in units of photons/s/cm^ x 10"^, and (panels 2-4) 
radio-band total flux density, fractional linear 
polarization (P%), and EVPA in the form of daily 
averages. Source-integrated data from the MOJAVE 
website are included for comparison. 

feet on the light curves. Simulated light curves from 
radiative transfer calculations based on this scenario 
may include multiple shocks within a single 'event' 
to explain the variability apparent in the radio band 
data. To reduce the number of free parameters, we 
make the simplifying assumption that multiple shocks 
within a single 'event' have the same obliquity. A re- 
cent refinement of our model has been the incorpo- 
ration of retarded-time effects. However, comparison 
of simulated light curves generated with and without 
retarded-time effects for test cases, shows that only 
small differences in the simulated flare shapes and 
spectral behavior result. 



4. Comparison of Data and Simulations: 
Method and Results 

Following the launch of Fermi in June 2008, a 
new phase of high-amplitude, centimeter-band activ- 
ity commenced. In this paper we study the 2009- 
2011 radio-band flare contemporaneous with 7-ray 
flaring. The combined radio and 7-ray data dur- 
ing this radio outburst are shown in Figure [H The 
weekly-binned photon fluxes shown in the bottom 
panel were obtained using ScienceTools-v9r27pl 
and P7SOURCE_V6 event selection. The LAT data 
were extracted within a 10° region of interest (ROI) 
centered upon the position of the target. We used 
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an unbinncd likelihood analysis (tool gtlikc) to deter- 
mine the photon fluxes by including in the model all 
of the sources within 15° of the target and by freez- 
ing the spectral index of all sources to the values in 
the 2FGL catalogue. For bins with TS< 10, we show 
a 2-a upper limit. A description of the UMRAO ob- 
serving and reduction procedures is given in [lT| . Sev- 
eral low level 7-ray flares were detected prior to the 
main flare which started near RJD 55170 (2009.93). 
Panels 2-4 show daily averages of the UMRAO total 
flux density and linear polarization monitoring obser- 
vations at 14.5, 8.0, and 4.8 GHz. Substructure is 
clearly apparent in the fractional LP light curve. The 
centimeter-band LP during the 7-ray flaring exhibits 
a systematic swing in the EVPA through 90°, and an 
associated increase in the fractional LP; this swing is 
particularly well-defined by the higher sampling rate 
at 14.5 GHz. The increase in fractional LP and or- 
dered swing in EVPA are the signature of the passage 
of a shock through the emitting region. 

A serious impediment to a detailed analysis of sin- 
gle dish light curves is the overlapping of multiple, 
subsequent flares, apparent as structure in the ob- 
served light curves. Guided by this structure and us- 
ing an outburst shape taken from a library of simu- 
lations as a template, we estimated the relative times 
and strengths of each shock to provide a first ap- 
proximation used as input for our detailed modeling. 
Our methodology for the deconvolution of the blended 
events differs from that adopted by other groups which 
comrnonly assumes a generic, exponential flare shape, 
e.g. Il2|. 

A simulation which reproduces the main features of 
the data is shown in Figure [31 These characteristics 
include the spectral evolution and amplitude of the to- 
tal flux density outburst, the magnitude of the change 
in fractional LP, and the frequency-dependent changes 
with time in the EVPA light curves. A comparison of 
the data and simulations is given in TablelH The prop- 
erties compared are: the range in the amplitude of the 
total flux density at the highest frequency (14.5 GHz), 
the evolution in the spectral index, 0:5, based on the 
total flux densities at 14.5 and 4.8 GHz, the ampli- 
tude of the change (minimum to maximum) of the 
fractional linear polarization, and the spectral change 
in the EVPAs during the evolution of the outburst. 
The latter is complex. As the flare starts, the EVPAs 
at 4.8 GHz are shifted from those at 14.5 and 8 GHz 
which track each other, but the EVPAs converge for 
the three frequencies just beyond the outburst peak. 
There is a null in the observed P% and an associated 
sharp variation in EVPA preceding RJD 55300 which 
is not reproduced by the model, however. In gen- 
erating these simulated light curves we incorporated 
three forward moving shocks oriented transversely to 
the flow axis, and to characterize the flow we adopted 
an optically-thin spectral index Q!i/ii„=0.25, a fiducial 
"thermal" Lorentz factor (associated with the random 
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Figure 3: Simulation of the radio band light curves 
during the 2009-2010 event. The computations have been 
carried out at 3 harmonically-related (\/3) frequencies 
which correspond to the UMRAO frequencies of 14.5, 
8.0, and 4.8 GHz; the symbols follow the convention used 
in Figures 1 & 2. The abscissa shows fractional time 
during the event modeled. 



motions of the emitting particles) 7c=1000, a low en- 
ergy cutoff Lorentz factor 7/ =50, a bulk Lorentz factor 
of the flow 7j=5.0, and an observer's viewing angle of 
9obs= 2°. See [l3| for definitions and explanations of 
these parameters. We tabulate in Table HIl the param- 
eters specifying the three shocks. The first shock is 
weak. The second shock is the main, central one, and 
it starts at around t=2009.7. The third shock com- 
mences around 2010.0, and it is this shock which coin- 
cides temporally with the brightest 7-ray flare. This 
shock is only 35% of the length of the second shock 
and, hence, a new jet component corresponding to 
this shock might not be expected to be resolved in se- 
quences of VLBI images. However, a new component 
is apparent in a MO JAVE imaging map obtained near 
this epoch as discussed earlier. Based on the Lorentz 
factor of the shock of about 12, and a viewing angle of 
2° we derive a Papp=8 c which is consistent with the 
MOJAVE results on component speeds. The shocked 
flow has a Lorentz factor of about 8. The model fit- 
ting also allows us to constrain the particle energetics; 
in a second source recently modeled, OJ 287, a lower 
energy cutoff of the 'thermal' Lorentz factor was re- 
quired to obtain a fit to the data. 



5. Discussion 

Multifrequency, single dish polarimetry data, com- 
bined with modeling, have provided sensitive probes 
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Table I Comparison of the Observed and Simulated 
Radio-Band Light Curves 



Property 


Observed 


Simulated 


Range in S(14.5) 


4-8 


3- 8 


as evolution 


flat, inv., flat 


flat, inv., flat 


Range in LP% 


1-3 


0-2 


EVPA behavior 


complex 


complex 



Table II Shock parameters 



Shock 


1 


2 


3 


obliquity (77) 


90 


90 


90 


sense 


F 


F 


F 


length (1) 


10 


11 


4 


compression (/t) 


0.8 


0.6 


0.6 



of the jet conditions associated with 7-ray flaring in 

this blazar. For fast flows seen at orientations of only 
a few degrees, as in this source, small changes in view- 
ing angle make significant differences in the simulated 
properties, most significantly in the amplitude of the 
fractional linear polarization. As a consequence, our 
LP data are a powerful constraint on this important 
parameter and an independent means for its deter- 
mination: viewing angle determinations obtained by 
comparing positions of core and jet features on VLBI 
maps are sensitive to the segment of the jet domi- 
nating the emission at the frequency used for the ob- 
servations. We note that the goal of this stage of our 
work has been to validate our assumed scenario, which 
for the source discussed here incorporates transverse 
shocks and rectilinear fiow. While we wish to explore 
other fits to this source and to model the variability in 
other sources in future work, real flows are complex, 
hydrodynamical entities. Jet curvature in particular 
on a variety of scales can affect the observer's inter- 
pretation of the observed variability. We caution that 
the simulations are highly sensitive to the flow de- 
tail, so, in reality, it may not be possible to improve 
dramatically upon the already-good fit described here 
for 0420-014. The method has yielded results consis- 
tent with VLBI determinations while provided more 
detailed information about the flow conditions where 
the flaring is produced. 
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